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Interferometry of non-Abelian edge excitations is a useful tool in topological quantum computing. 
In this paper we present a theory of a non-Abelian edge state interferometry in a 3D topological 
insulator brought in proximity to an s-wave superconductor. The non-Abelian edge excitations 
in this system have the same statistics as in the previously studied 5/2 fractional quantum Hall 
(FQH) effect and chiral p-wave superconductors. There are however crucial differences between the 
setup we consider and these systems, like the need for a converter between charged and neutral 
excitations and the neutrality of the non-Abelian excitations. These differences manifest themselves 
in a temperature scaling exponent of —7/4 for the conductance instead of —3/2 as in the 5/2 FQH 
effect. 



PACS numbers: 03.67.Lx 71.10.Pm 73.23.-b 74.45. +c 

I. INTRODUCTION 

One of the most promising tools in topologi- 
cal quantum computing-^ is non-Abelian edge state 
interferometryJSEEl Its main idea is that moving a frac- 
tional excitation (anyon) existing at an edge of a topolog- 
ical medium around localized anyons in the bulk allows 
to extract information about the state of the latter. The 
theory of edge state interferometry was initially devel- 
oped for Ising anyons in the 5/2 fractional quantum Hall 
(FQH) state and p-wave superconductors,'^'^' building 
on earlier work on FQH systems!^ Recent experiments,^ 
which provide evidence for non-Abelian braiding statis- 
tics in the 5/2 FQH state (see the detailed discussion in 
Ref. [TUJ) are using this method, and it is generally con- 
sidered the most promising way to measure the state of 
topological qubits. 

We present a theory of non-Abelian edge state in- 
terferometry of the Majorana modes existing at the 
surface of a 3D topological insulator brought in con- 
tact with an s-wave superconductor and a ferromag- 
netic insulator.^ The main difference of an interferome- 
try setup in this system, as compared with 5/2 FQH in- 
terferometer, is t he nee d for an additional "Dirac to Ma- 
jorana converter" GHH] This element is required because 
unlike in the FQH effect the edge excitations near a su- 
perconductor carry no charge and thus allow no electric 
readout. This converter initially transforms the charged 
excitations injected from a current source into superposi- 
tions of two neutral excitations existing at different edges 
of the superconductor. Later another converter recom- 
bines a pair of neutral excitations exiting the interferom- 
eter into a charged particle, either an electron or a hole, 
that can be measured as a current pulse. The difference 
between the two systems is summarized in Fig. [I] The 
"Dirac to Majorana converter" is not available in chiral 
p-wave superconductors, since the chirality of the neutral 
edge modes is then set by time-reversal symmetry break- 
ing in the condensate, and not by the external region of 
the system (magnet). Such a limitation combined with 



the absence of charged modes makes electric readout of 
interferometry experiment much less viable in a chiral 
p-wave superconductor. 

The description of the "Dirac to Majorana converter" 
using single particle formalism was done in Refs. 1121131 
The qualitative description of the non-Abelian Fabry- 
Perot interferometer was presented in Ref. [T31 In this 
paper we use conformal field theory (CFT) to describe 
and analyze the non-Abelian excitations following Ref. 6 

An important difference between the systems is the 
following: In the 5/2 FQH effect the charge density and 
accordingly charge current of anyons may be defined lo- 
cally, since anyons have charge e/4 or e/2 in this system. 
Excitation of charge e* has an energy cost of e*V for 
being created in the system. This energy cost provides 
a natural cutoff for the current, whereas in the super- 
conducting systems due to the absence of charge in the 
edge excitations the only cutoff is set by the finite tem- 
perature. The neutrality of the edge excitations does not 
only mean that a finite voltage does not provide a cut- 
off for the conductance, but also results in a different 
temperature scaling exponent of the conductance. In the 
topological insulator setup the conductance diverges at 
low temperatures as T ~ T~ 7 / 4 , while in the FQH setup 
it goes as T - T~ 3 / 2 . 

The experimental requirements for a realization of edge 
state interferometry in topological insulators were dis- 
cussed in Refs. 12 and 13. An additional requirement 
for non-Abelian interferometry is the need for a suffi- 
ciently high amplitude of the vortex tunneling, \ a ~ 
exp(— yjEc/Ej), with Ej the Josephson energy and Eq 
the charging energy. It is non-negligible only if the su- 
perconducting islands in the system have small capacitive 
energy £?c- 

The outline of this paper is as follows: In Sec. [IT| we 
introduce the effective model that we use to describe the 
fcrmions that propagate along magnetic domain walls 
and the superconducting-magnet domain walls. In par- 
ticular we introduce the representation of these fcrmions 
in terms of Majorana fields, which we use later. In Sec. |III| 
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FIG. 1: Edge state Fabry-Perot interferometer in the 5/2 
FQH system (top panel) and in a topological insulator/s-wave 
superconductor heterostructure (bottom panel). The charge 
is transferred locally at the tunneling point in FQH effect, 
and is only well-defined in the ferromagnetic domain walls 
(i.e. the leads) in the topological insulator setup. Regions 
labeled S, Mf , and Mi denote parts of topological insulator in 
proximity of a superconductor and of ferromagnetic insulators 
with different polarizations. Grey circles in the middle of the 
central island are Majorana bound states forming a Majorana 
qubit, which can be measured by the interferometer. 



we review the linear response formula that we use to cal- 
culate the non-local conductance, the experimentally rel- 
evant quantity that we are interested in. In Sec. |IV| we 
give a detailed account of the perturbative calculation 
of the conductance, and we consider the most interesting 
case of vortex tunneling in Sec.[V] In Sec. |VI| we show how 
the proposed setup can be used to measure the fermion 
parity (and hence the topological charge) of the Majo- 
rana qubit that is stored in a pair of bulk vortices. Our 
conclusions are to be found in Sec. |VII| We provide a 
detailed description of the formalism that we use to de- 
scribe the peculiar vortex field in the appendices. 



II. CHIRAL FERMIONS 

A. Domain wall fermions 

It is known that there exists a single chiral fermion 
mode on each mass domain wall in the 2D Dirac equa- 
tion. This mode is localized near the domain wall but is 
allowed to propagate along the domain wall in only one 
direction (hence the name chiral). This is most easily 



seen using an index theorem that relates the difference in 
a topological number (N 3 in the language of Ref. [15]) be- 
tween the two domains and the difference in the number 
of right- and left-moving states that live in the domain 
wall. 15 In the ferromagnetic domain wall that we are in- 
terested in the change in N3 across the domain wall is 
±1. If the domain wall is also abrupt enough then only 
one chiral fermion exists in the domain wall. 

A similar argument can be made using the Dirac- 
Bogoliubov-de Gennes (BdG) equation with gaps gen- 
erated by the superconducting order parameter A. In 
the case that we consider (s-wave pairing) A3 is zero if 
the gap is dominated by the superconducting gap |A| and 
non-zero (±1) when the gap is of ferromagnetic charac- 
ter. Because of the double counting of states in the BdG 
equation this implies that | of a chiral fermion state ex- 
ists on a superconducting-magnetic domain wall. This is 
exactly the number of degrees of freedom that is encoded 
in a chiral Majorana fermion field. 

Alternatively one can argue for the existence of these 
states by solving the BdG equation explicitly for certain 
simple domain wall profiles or use k • p theory.^ We now 
proceed to a theoretical description of these states. In 
particular we will see that it is fruitful to describe both 
kinds of domain walls in terms of Majorana fields. 



B. Theoretical description 

In the leads (ferromagnetic domain walls), where the 
superconducting order parameter vanishes, the system 
consists of a single normal edge state which propagates 
in only one direction, i.e. a single chiral charged mode. 
This can be described by a complex fermionic field ^(x) 
with Hamiltonian 



H(t) 



1 

2^ 



dx :&(x)[vp x - p(x, t)]$>(x) : 



(1) 



Here : : denotes normal ordering. We use units such 
that Ti — 1 unless specified otherwise. The kinetic energy 
operator vp x is defined as 



.d x v(x) - v(x)d x 



~i^v(x)d x ^v(x), (2) 



where we have introduced the spatially varying velocity 
v(x) in a symmetric way such that vp x is a Hermitean 
operator. The stationary (energy E) solution to the time- 
dependent Schrodinger equation corresponding to Eq. ([IJ 
for zero chemical potential p, — is 



v(0) 



v(x) 

This implies that 



exp ( iE 



dx' 



(0,0)) = 



v ( x ') 

v(zMo)]- 1 / 2 



*s(0,0). (3) 



a + i[t - J* dx> /v(x')] 



(4) 
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where a is a short time cutoff which should be taken to 
zero, ff the velocity v is constant the result simplifies to 



u(*(x,t)* T (0,0)) = 



f 



f 



a + i(t — x/v) a + 



(5) 



The normalization in Eq. (JT|) is chosen to yield this result 
without any extra normalization factors. Note that it 
implies (in the limit a — > + ) that the anti-commutation 
relation for the field is {#(x), &(x')} = 2nd(x - x'). 

An important consequence of the chiral nature of the 
excitations is that the correlation functions only depend 
on the difference of the Lorentz time u = t—x/v. Accord- 
ing to Eq. Q the same is true also for a spatially varying 
velocity with the proper interpretation of the length dif- 
ference. Because of this property we will mostly work 
with a spatially homogeneous velocity that we will set 
to unity (v = 1) in the following calculations. It is also 
useful to go from the Hamiltonian to the corresponding 
Lagrangian 

L = ^ j ' dx:¥{x)[id t -vp x + n{x,t)]V(x):, (6) 

since the coupling to the gauge field is most transparent 
in this formalism. 



C. Majorana fermion representation 

We can decompose ^>(x) into two independent Majo- 
rana fields ip(x) = t/)'(x) and ip'>(x) = ip'( x ) as 



%A[x,t) 



-[iP{x,t) + iiP'{x,t)]. 



(7) 



The anti-commutation relations of the Majorana fields 
are {ip{x),i/j{x')} = {tp'(x), f{x')} = 2ir5{x - x'), and 
{tp^x), ip'(x')} — 0. In terms of tp and ip' the Lagrangian 
becomes 



1 



L=—J dx[:ip(x)(id t - vp x )ip(x) : 

+ :ip'(x)(id t - vp x )ip'(x) :] 

+ — / dxF(x,t)v(x)ij'(x)ip(x), (8) 

where F(x,t) depends on the phase A(x,t), i.e. it is 
gauge dependent: 

- eF(x, t) = - -^d t A(x, t) - d x A(x, t) . (9) 

Note that this means that a time-independent spatially 
varying chemical potential can be gauged away up to 
possible boundary terms. 

One of the most interesting features of the system that 
we consider is that the two Majorana fields that appear 
in this action can becomes spatially separated when a su- 
perconducting region is sandwiched in between the two 



magnetic regions in a magnetic domain wall as discussed 
previously. Thus the action in Eq. ^ can be used to 
describe the setup in Fig. [2] in which the two Majorana 
fields ip and are spatially separated inside of the in- 
terferometer. It is important to remember that the co- 
ordinate systems of the two fields are different in this 
representation. 




FIG. 2: Free fermion propagation setup. The two Majorana 
modes ip and ip' are spatially separated by the superconduct- 
ing region. Thus the effective propagation length from in to 
out can be different for the two modes, i.e. L' ^ L. 

From the Lagrangian and the coupling to the gauge 
field we now identify the charge current operator as 



J(x) 



:$Ux)%(x):= ^v(x)ip'(x)xp(x). (10) 

Z7T Z7T 



This form of the current operator in terms of the Ma- 
jorana fields is very important for the following calcula- 
tions. It is only well-defined if the two Majorana modes 
are at the same position in space, hence there is no cou- 
pling to the electric field inside of the interferometer 
where the two Majorana wires are spatially separated. 
This is also an important difference between the FQH 
setup where local charge current operators can be de- 
fined at the tunneling point contacts. This simplifies the 
calculation because the local charge transfer is directly 
related to the measurements done far away. In our sys- 
tem we don't have this luxury and must consider the 
leads explicitly. 



III. LINEAR RESPONSE FORMALISM FOR 
THE CONDUCTANCE 

If we write the Lagrangian in Eq. pj as L = L — H'(t), 
where the term on the last line is 



H'(t) = dxJ(x,t)F(x,t), 



(11) 



we are in the position to use the standard linear response 
Kubo formula,^ to calculate the conductance tensor T. 
Following Ref. [17] we introduce an AC chemical potential 
localized in the source lead, which we take to have coor- 
dinates x < 0. We choose a constant gauge A(x, t) = A 
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so that F(x,t) = -e(-x)cos(nt)e- s ^V/v(x)^ The 
conductance T is defined as the magnitude of the in- 
phase current divided by the applied voltage difference 
V. Following the usual steps, with the current operator 
in Eq. ( 10 1 and assuming that the two Majorana modes 



the expressions for the Green's functions into Eq. ( 14 1 we 
obtain 



are independent, we obtain the formula 



e 

■nh fi 



/>oo POO 

lim / dt' / ^Im[G> l G > , i ,]cos(nt)e" 



til 



(12) 



Here we have reintroduced the correct units of conduc- 
tance e 2 /h. We have also used the fact that in a chiral 
system the response in the region x > to a spatially uni- 
form extended source x < at a particular time t' = is 
equivalent to the response to a point source at x = that 
is on for t' > 0. The important quantities to calculate 
are the Green's functions 



G j'i' 



(lKl/,f)V(0,O)=<V'ilfc>> (13a) 
^'(y',t)^(0',t'))^(^,). (13b) 



Here the indexes i and j are shorthands for the coordi- 
nates of the source (0,i') and current measurement (y,t). 
Similarly for the primed coordinate system, which is typ- 
ically not the same in the setups that we consider as dis- 
cussed previously. 

Because the correlation functions only depends on t— t' 
it is possible to perform the integral over t + t' in Eq. ( 12 ) 
explicitly, the resulting expression is 



e 2 f°° 

r = -- h J o dtiMG>G> t ,]t, 



(14) 



where it is understood that the source term is taken at 
t' = 0. Here we have also used the fact that the correct 
limit is to take 5 — > + first and then O — > 0. Because we 
are interested in the finite temperature result the cut-off 
provided by the thermal length is enough to render the 
expression convergent. This is the master formula that 
we will use to calculate the conductance in the following. 



A. Free fermion propagation 

If both Majorana modes propagates freely (the setup 
is sketched in Fig. |2| we can use the finite temperature 
propagator 



ttT(L - L') 



h sinh[7rT( J L - U) 



(16) 



in the limit a — > + . This formula agrees with the lin- 
ear response limit of the the result obtained with the 
scattering formalism in Ref. 12, and shows how the path 
difference enters in the finite temperature case. 

To obtain the response in the source lead we take the 
limit L' — > L with the result that 



r = 



(17) 



This is the expected (and correct) result for a system 
with one propagating channel. If L' =/= L we also obtain 
Eq. (17 1 as long as T\L — L'\ <C 1, in the zero temperature 
limit the result is thus independent of the path length dif- 
ference. The E q. ^Tf\ agrees with the limit V — > of the 
previous resultsp^ 1 ^ which were based on the scattering 
formalism. 



IV. PERTURB ATIVE FORMULATION 



In tunneling problems we want to calculate the Green's 
function Gj* = (tpjipi), where ipi and ipj live on different 
edges of the sample, in the presence of a perturbation 
5 H that couples the two edges. Assuming that the sys- 
tem is in a known state at time to, we may express the 
expectation value in the interaction picture as 



G> = (tf(i O) i)^(t)tfMM(0M0,t )}- 



(18) 



Here U(t, t') is the time evolution operator in the inter- 
action picture. For t > t' it is given by the familiar time- 
ordered exponential U(t,t') = Texp[— i 

In the following we will assume that the average at 
t = to is a thermal one at temperature T. A perturbative 
expansion is obtained by expanding the time-ordered and 
anti-time-ordered exponentials in this expression in pow- 
ers of 5H . This procedure is equivalent to the Schwinger- 
Keldysh formalism, which in addition provides a scheme 
to keep track of whether one is propagating forward or 
backward in time. We will also assume that the per- 
turbation was turned on in the infinite past, i.e. we set 
to = — oo. 



c > i= 1 _ *T 
31 Zj'i' simrT[a + iiij'i'] 

7rT 

= n5(uj<i>) -iV , r, (15) 

a^o+ smh(7rTwjv) 

where Uj'i' = t — L' . The Green's function of the other 
edge is given by the same expression with L (the ef- 
fective length of propagation) instead of L' . Substituting 



A. Fermion tunneling 

As a warm-up for the vortex tunneling calculation we 
will now consider the simpler case of fermion tunnel- 
ing, which we describe by a tunneling term H^(ti) = 
1X^2^1 /(27r)P Here ipi (^2) is located at the tunneling 
point at the upper (lower) edge. The system and the co- 
ordinate convention we use are sketched in Fig. [3] The 



5 



leading contribution to conductance comes at first order 
in the tunneling amplitude A^. After a straightforward 
expansion and collection of terms we obtain 



Aw, 



2tt 



dti{^,^i}<^> + 0(A4). (19) 



Here we have used the fact that the two groups of 
fermions on different edges, i.e (ipj,ip2) and {ipi,ip\), 
are independent. It is straightforward to evaluate this 
expression using Eq. ( |15[ ) together with ipi} — 
2ttS(uu), and {ipj, 1P2} = ^5{u 2 j), where uu = ti — L iop 
and U2j — t\ — t + £bottom- Because of the geometry 
of the problem the second term on the right hand side 
of Eq. (19) vanishes due to causality (the Lorentz time 
arguments never coincide). The Green's function to 
leading order in tunneling strength is therefore 



VORTEX TUNNELING 



The main focus of this paper is to study how the 
tunneling of a pair of vortices can effectively transfer a 
fermion, and hence give a contribution to the conduc- 
tance. Schematically the vortex tunneling term can be 
written as 



(22) 



where the index t (b) denotes the top (bottom) edge. As 
it stands this term is not well-defined without more in- 
formation about the two spin fields a, this is discussed in 
great detail in Ref. H) We provide a detailed description 
of the formalism that we use to deal with this issue in 
the appendixes. 



simrT[a + i(t - L)]' 



(20) 



Coordinate conventions 



where L — L top + ^bottom is the effective propagation 
length of the Majorana fermion. Using the result of 
Sec. |III A| we then find that the conductance of this setup 
is 



\ e 2 

A,/, , 



(21) 



at T = 0. Once again this result agrees with the zero 
frequency, zero voltage limit of the results obtained with 
the scattering method in previous workP^EH 
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FIG. 3: Top panel: fermion tunneling setup. The coordinate 
conventions used in Sec. |IV A] are shown in the bottom panel. 



To have a well-defined prescription for the commuta- 
tion relation of fields on different edges we will treat the 
two edges as spatially separated parts of the same edge. 
This reasoning has been employed in a number of works 
studying tunneling in the FQH effect, see for example 
Refs.lT9landl2"Ul This approach leaves a gauge ambiguity: 
should we choose the bottom edge to have spatial coor- 
dinates smaller or larger than that of the top edge? The 
correct choice is fixed by noting that the current operator 
at the source should commute with the vortex tunneling 
term at equal times because of the locality and gauge 
invariance. A similar argument can be made considering 
the current operator at the measurement position before 
the information about the tunneling event has had time 
to reach it. Since we want the vortex tunneling event to 
commute with fermions on the reference edge at all times 
we are forced to use the coordinate convention shown in 
Fig. [4] in which the spatial coordinates on bottom edge 
are always larger than those on top edgepf The vortex 
tunneling then corresponds to changing the phase of the 
superconducting order parameter by ±27r to the right of 
the tunneling point in the figure. 

In addition it is convenient to introduce an even 
more compact notation. We denote ip t (0,—L t ) = ip^, 
tpb(t, AL + L b ) = ipj, a t (ti,-x t ) = 01, a b (h,AL + Xb) = 
02, Ct(t2,xt) = 03, 06 (i2, AL— Xb) = 04- The two tunnel- 
ing terms in the Hamiltonian are then written as \ a T\2 
and A CT T34. The modification needed to allow for dif- 
ferent tunneling amplitudes A ct l and A a ^ at the left and 
right tunneling points (see Fig.[T| is straightforward. The 
"Lorentz times" u for right- movers are u = t — x. We 
use additional short-hand notations u a p = u a — up and 
s a fj = sign(u Q — up). The Lorentz times of the six oper- 
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FIG. 4: Top panel: independent coordinate system for the 
two edges. Bottom panel: coordinate system in which the 
two edges are treated as spatially separated parts of the same 
edge. This allows us to correctly capture the commutation 
relations of the fields on different edges in the relevant limit 
AL -> oo. 



ators used in the calculation are 



"2 



L L 

t\ + x t 

h - x t 

t - L b - AL 

ti — Xb — AL 

h + x b - AL. 



(23) 



Taking the limit of large spatial separation AL — » +oo 
we see that s,y = 1. Accordingly, in this limit also Ski = 1 
for any k 6 {i, 1, 3} and I £ {j, 2, 4}. 

In the following perturbative treatment we will assume 
that ti >t\. This means that to calculate the full Green's 
function G^ we should sum over the four processes for 
which the first and the second vortex tunneling events 
happen at the right or the left tunneling point. The am- 
plitudes of the two processes in which vortex tunneling 
events occur at different points are related by changing 
Xt — > — Xt and Xb — > —Xb- Likewise the amplitudes of the 
processes in which both evens occur at the same tunnel- 
ing point can be obtained from the amplitude of the pro- 
cess with vortex tunneling at different points by setting 
Xt = Xb = and setting L t — > L t ± xt and Lb — > Lb ± Xb- 



B. Perturbative calculation of G 

In the appendices we demonstrate how one can eval- 
uate the averages of the contributions to the integrands 
generated in the perturbative expansion of G^. The tech- 
nically simplest way of performing the calculation is to 
use the commutation relation between fcrmions and tun- 



neling terms [see Eq. (B8|] 

^12^3 = Si 3 S2 3 1p 3 T 1 - 



to transform the correlation functions into one of the two 
forms in Eq. 
AL 



(B9I 



The limit of large spatial separation 
oo can then be taken using Eq. (BIO I. Finally we 



use the functional form of the correlation function of a ijj 
and two cr's that is fixed by conformal invarianceP^ 



3/8 



V2zl 



1/2 1/2' 

Z 3i 



(25) 



The result of this calculation is the same as the limit 
AL — > 00 of the full six-point function that can also be 
calculated using bosonization and a doubling trick, see 
App.0 

The first non-vanishing contribution to the fermion 
propagator G > comes at second order in the vortex tun- 
neling term. It is then convenient to divide the interme- 
diate time integrals into different regions. We will use 
the following labeling conventions: (a) t\ < t% < 0, (b) 
t\ < < t 2 < t, and (c) < t% < t 2 < t. We now calcu- 
late the contribution to the integrand from each region 
separately. 

Let us first consider the interval t\ < t 2 < 0. By 
straightforward expansion, and using the exchange alge- 
bra we obtain the integrand in this region 



I{ a ) = (V^TkTb) + {T^T^ipj^) 

- (T 3 4,ip 3 ipiTi2) - (T 12 ipjtpiTsi) 

= SiiSi2(Si 3 Si4 - S3jS4j){l/}jT3iTi2lpi) 

- SijS2j(Si3Sii - S 3 jS4j)(^jT 12 T 3 4lpi 



(26) 



The minus signs are generated when the two tunneling 
terms are on different Keldysh branches, i.e. when one 
comes from evolving forward in time and one backwards. 
We can simplify this expression further by noting that 
because of the geometry we always have Si 3 = Sn = 1 in 
this region. Thus 

1(a) = ^ = (1 + S]4)((lpjT 3 4T 12 ^i) + Sj2{i>jTi2T34lpi)). 

(27) 

Let us now consider the interval t\ < < t 2 < t. We 
denote the contribution to the integrand in this region 
by I n,y Expanding we get 

1(b) = (TizT^ipjipi) + (ipjTMipiTiz) 
- (TiztpjTuipi) ~ (T 34: ip 3 ipiTi2) 
= ...=/>. (28) 

To see that we get the same expression as in region (a) we 
have used the fact that Sn = 1 in this region. Performing 
the same calculation as in regions (a) and (b) for the 
interval < t\ < t 2 < t we find that also in this region 



1(c) - I > , 



(29) 



(24) 



and hence we can use I y throughout all regions. Us- 
ing cluster decomposition (i.e. taking the limit of spatial 
separation) and the explicit correlation functions we get 
the expression for the integrand. Putting back the inte- 
grals and the strength of the tunneling term we obtain 
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the leading term in the perturbative expansion of the 
Green's function 



G ; 



3_ 

2 3 / 2 



dti \ dt 2 



(1 + Sji ) 



(\z j2 \\z j4 \y/2(z 3iZli )i/2 



[(1 + s j2 ) Re(zU 8 z% 8 ) - (1 - s j2 ) Im(4 /8 z 4 3 2 /8 )]. (30) 

Note that this expression is a short form that includes 
a sum of many terms, it is valid for real times only and 
the analytic structure of the the Green's function is not 
apparent. It is useful to shift the time-coordinates tj = 



t 



G> = 



Xb - 



for j = 1,2. The resulting expression is 

f Sl 1 
dsi I ds 2 



V^Jo """Jo (kj 2 |M) 1/2 (^:u) 1/2 
[(1 + s 32 ) Re(zU 8 z% 8 ) - (1 - s j2 ) Im(z 3 3 / 8 4 /8 )], (31) 

where 



u j2 

Uu 

U31 
U42 

i 



si 



2x b 
s 2 

t + x t -x b ~ 
t - x t - x b - 
si - s 2 - 2x t 
si — s 2 + 2x b 
t — Lt — L b . 



si 
s 2 



(32) 



Note that the dependence on the parameters t, L tl and L b 
only enters in the combination t. The analytic structure 
is much more transparent in this equation. For tunneling 
at the same point, i.e. x t = x b = 0, we always have 
Sj 2 = 1 and the result simplifies to 



A*V2cos(y 



dsi / ds 2 
Jo 

Ni| 3/4 



(33) 



From this expression we see that Re[G > ] ^ only for 
times such that t > L t + L b . Since Re[G > ] is proportional 
to the retarded Green function G R , this is a reflection of 
the causality of the theory: information has to have time 
to propagate through the system for G R to be non-zero. 

The Green's function G > has a singular part that is 
given by 



G >^ A 2 T -3/4 Hlog | e | + 7re(0]j 

with 0(x) the Heaviside step function and 

£ = T(t - L t - L b - x t - x b ) < 1. 

C. Conductance 



(34) 



(35) 



Substituting the propagator in Eq. (15 1 for the refer- 



ence edge into the expression for conductance in Eq. ( 14 ) 
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CN b 



5r 



2; 

oL 
0.0 




FIG. 



2.0 2.5 



Normalized conductance (h/e 2 )T /A„ 



T" 7/4 F[0, 0, {L t + L b )T, L'T] of a single quantum point con- 
tact due to vortex tunneling as a function of temperature. 
The parameters of the setup are L t — Lb = L' . 



we obtain 



2 /h 



dtV 



Tt 

sinh(7rT'Ujv) 



Re[G>]. (36) 



Together with Eq. (31 ) this expression provides a closed 



expression determining the contribution from each pro- 
cess to the conductance, which may be directly evaluated 
numerically. Since G > only has a logarithmic divergence, 
the short distance cutoff a may be directly set to zero in 
this expression. By substituting the singular part of G > 
into Eq. ( 36 ) one can see that the conductance contribu- 



tion is a continuous function of all the parameters of the 
problem. It may be written as 



F 



LB. 



\ 2 a F[x t T,x b T, {L t + L b )T, L'T] 

J>7/4 



(37) 



with F a universal continuous function. In the low tem- 
perature limit, when all of the arguments of F are small, 
the contributions to conductance from vortex tunneling 
at different points Tll, T rr, Tlr, and Trl are all equal 
to each other and to 



r = 



e 2 A^O, 0,0,0) 
h T 7 / 4 



(38) 



with F(0, 0,0,0) « 1.5. In the other limit, when either 
\x t + x b \T > 1 or \L t + L b - L'\T > 1 the function F 
is exponentially small, or in other words conductance is 
suppressed due to thermal averaging. We have evaluated 
the conductance of a single point contact due to vortex 
tunneling numerically with the result shown in Fig. [5] At 
low temperatures T x T 7 / 4 — > constant as expected, and 
at high temperatures T ~ exp(— T\L' — L t — L b \). 

The scaling exponent of conductance —7/4 is different 
from —3/2, the exponent of tunneling conductance in 
the 5/2 FQH effect. This naturally follows from the very 
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different mechanisms of conduction in the two systems: 
current is carried by charged modes in 5/2 FQH system, 
while "Dirac to Majorana converter" forms current in 
topological insulators. 



VI. QUASICLASSICAL APPROACH AND 
FERMION PARITY MEASUREMENT 

The most interesting application of the interferome- 
ter setup with vortex tunneling is that it allows for the 
detection of the fermion parity of th e sup erconducting 
island between the two point contacts! 3 * 4 * 5 * This is possi- 
ble because vortices acquire a phase of it when they are 
moved around an odd number of fermions P^l In the sim- 
plest case, when there are only two bulk vortices in the 
central region, as shown in Fig.[l] the interferometric sig- 
nal reads out the state of the qubit formed by the bulk 
vortices. 

Without loss of generality we consider the case of two 
bulk vortices that are situated in between the left and 
the right tunneling regions. From the point of view of 
the electronic excitations the bulk vortices can be de- 
scribed by two localized Majorana bound states,^ with 
corresponding operators j a and 7&. To describe the ac- 
tion of the vortex tunneling term on these excitations 
we include, following Ref. 1241 an extra term P a b = i^/ a jb 
in the left tunneling operator. This operator captures 
the property that upon changing the phase of the order 
parameter in the superconductor by ±2n the Majorana 
modes localized in the vortex cores gains a minus sign. 

In the absence of bulk-edge coupling the fermion par- 
ity of the vortex pair is a good quantum number that 
does not change with time. In that case the extra term 
that is added to the left tunneling term P a i, measures the 
fermion parity of the qubit defined by 7 a and 7^. This 
means that we can replace P a b (— l) n/ , where n/ is 
the number of fermions in the two vortices. In the second 
order calculation this factor enters only in the contribu- 
tions where one vortex tunnels at the left tunneling point 
and one at the right, so the total conductance is equal to 



r = r LL + r RR + (-iy i f (r LR + r RL ). 



(39) 



The expressions for the Vs were calculated in the previ- 
ous section. The effect of bulk-edge coupling is presum- 
ably similar to the case of the 5/2 FQH effec t that has 
been studied in great detail recentlyp3BSB5EZl 

The phenomenological picture of the non-Abelian in- 
terferometry presented in Ref. [T3] can be summarized in 
the following way. First an incoming electron is split 
into two Majorana fermions when it approaches the su- 
perconductor. Next one of these Majorana fermions is 
further split into two edge vortices, or a excitations. The 
edge vortices tunnel at either of the two point contacts, 
and recombine into a Majorana fermion again. Finally 
two Majorana fermions combine into electron or a hole as 
they leave the superconductor. At zero voltage any dy- 
namic phases are prohibited by electron-hole symmetry, 



so the outgoing current may be written as 

I = jV [\ 2 aL + \ 2 aR + 2{-l) n f\ aL ~\ aR ] , (40) 

where \ aa (with a = L, R) is an effective vortex tunneling 
amplitude (here we allow for different vortex tunneling 
amplitudes at the left and right tunneling points). 

Comparing Eqs. (l38])-((39l with Eq. (|4Q| we see that 



at low temperatures the effective vortex tunneling am- 
plitude is equal to 



Ka = XaaT- 7/8 ^F(0,0,0,0). 



(41) 



Once this identification is done, the quasiclassical picture 
is directly applicable given that 1/T is much larger than 
the characteristic length of the system and the second 
order perturbation theory still holds (X aa <!)■ 



VII. CONCLUSIONS 

In this paper we have introduced a theory for a non- 
Abelian interferometer on the surface state of a 3D topo- 
logical insulator brought in proximity to an s-wave super- 
conductor. This theory uses CFT to describe the vortex 
field following Ref. [6] and is an extension of the earlier 
qualitative discussion in Ref. [131 In particular we showed 
that if the temperature is low and tunneling is sufficiently 
weak, it is possible to introduce an effective tunneling am- 
plitude of vortices according to Eq. (41 1. This justifies 



the simple quasiclassical description of vortex tunneling 
used in Ref. |T3] 

Because the vortex tunneling term is a relevant op- 
erator, the perturbative treatment is only valid at high 
enough temperatures. This statement is reflected in the 
divergence of conductance F ~ T~ 7 / 4 . The scaling ex- 
ponent —7/4 is different from the tunneling conductance 
scaling exponent —3/2 of the 5/2 FQH setup in the linear 
response regime due to the different structure of current 
operators in the two systems. 
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APPENDIX A: VORTEX TUNNELING TERM 

In this appendix we show how one can calculate the 
amplitude for transferring a fermion between the two 
edges in terms of two vortex tunneling events using 
bosonization with the help of a doubling trick. This 
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is an old technique that goes back to the seventies,^ 
which is now textbook material! 22 * 29 ^ In the appendices 
we use the condensed coordinate conventions introduced 
in Sec. |V A| but we'll keep the gauge choice implied by 
the sign of Sij unspecified. 



1. Non-chiral extension of the system 

The logic of the procedure can be motivated as fol- 
lows (see also the construction in Ref. l2"¥)l . We are in- 
terested in the tunneling of a chiral Majorana fermion 
between two edges of a sample (cf. Fig. [3]). Because of 
the fermion doubling feature it is convenient to enlarge 
the system by adding an additional counter-propagating 
chiral Majorana fermion. These two copies can then be 
described as the continuum limit of a lattice model of 
local Majorana fermions (described by lattice operators 
7/ = 7;) that are allowed to hop to their nearest neigh- 
bors: 



H 



2N 

-t^ilili+i- 

1=1 



(Al) 



2N 



t/4 



11- 



The fermion parity operator is then P = \\ 
This system is known to map onto the (quantum) Ising 
chain in a transverse field at criticality (see e.g. Ref. l29]l . 
which is also equivalent to the classical 2D Ising model 
at its critical point. In the Ising model there are spin and 
disorder fields that are non-local in terms of the lattice 
fermions. It is easy to write down explicit expressions 
for the spin and disorder operators in terms of a string 
of Majorana fermions on the lattice, for example 



2.7 



0"2i+l°"2j + l 



&2i+lL>-2j 



n eW4 ^ 

l=2i+l 
2J-1 

n 

l=2i 

2j-l 



(A2a) 
(A2b) 
(A2c) 



It is clear from these expressions that a cr/i term changes 
the fermion parity of the system whereas aa and [i\x do 
not. 

Now we are not interested in the lattice theory itself 
but rather the low-energy theory which is obtained in the 
continuum limit of the lattice model. This limit is known 
to map onto the Ising CFT. This is a thoroughly studied 
system and we can hence rely on results from the large 
literature on this topic. 

In particular, on the lattice we know that a vortex 
tunneling term has to be of the form <J\<72 or H1H2, oth- 
erwise the fermion parity is changed. Furthermore, from 



the operator product expansion of the Ising CFTP^EHl 



0102 



7 — m + 7:(zi2Zi2) 3/8 i^2, (A3a) 

(zi2Zi 2 ) 1/8 2 

~ 7 _ s 1/8 - ^(^i22i 2 ) 3/8 iV'2V'2, (A3b) 

(2l22:i2) 1/S 2 



we see that a pair of cr's (or a pair of /z's) can change the 
parity of right-movers. Since our tunneling term is not 
allowed to do this we take the tunneling term in the non- 
chiral system to be T12 oc criO"2+MiA t 2- Clearly the parity- 
changing term is canceled with this choice. Another way 
of putting this is to say that this combination enforces 
the tunneling term to be in the identity channel. 

It is known that two independent copies of the Ising 
model can be bosonized using Abelian bosonization I22I2S1 
It is then a straightforward calculation (using for example 
the explicit expressions in the appendix of Ref. |3"CT|) to 
show that the doubled tunneling term can be bosonized 
as 



T\2T'\2 



(4> 1 ~4> 2 \ (<i>x-<h\ 
\—^) cos {—^)- 



(A4) 



It is important to note that the primed system is an inde- 
pendent copy of the system in this expression, and that it 
is introduced as a trick to allow for a simple calculation 
of various correlation functions. 



2. Prom non-chiral back to chiral 

Since we are only interested in the right-moving part 
of the tunneling term we would like to get rid of the 
left-moving part in the last equation. Because of the 
factorization of the right- and left-moving parts we are 
allowed to use 



T 12 T' = cos 



(91 ~ <P2\ 



(A5) 



as the doubled tunneling term in the chiral system. 
Here the cosines are to be understood as shorthands for 
cos(a - b) = (e la e- ib + e - ia e lb )/2. The exponentials in 
these expressions are actually dimensionful vertex opera- 
tors, see e.g. Ref. I3T1 for a detailed discussion. With this 
representation together with the bosonized representa- 
tion of the Majorana fermion in the unprimed system 



ipi = \/2cos(<fo), 



(A6) 



and the standard bosonization formula (which holds if 



Si=i a i — 0, otherwise the expectation value vanishes) 



n 



l<i<j<N 



with 



sm[7rT(a + iui 
7rT 



(A7) 



(A8) 
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we can in principle calculate any correlation function us- 
ing the bosonization formalism. In particular we can cal- 
culate the full six-point function including two ?/>'s and 
two tunneling terms. This will be done in the next sub- 
section, but let us first check that the representation re- 
produces known results for the 2-, 3- and 4-point func- 
tions. 

Let us first consider the vortex 2-point function. This 
is calculated via 



(T 12 ) 2 = (T 12 T[ 2 ) = 



1 

~T/4' 

h 2 



(A9) 



Taking the square root we obtain the correct result for a 
field with dimension 



16 



(Tl2) = 



1/8' 



(A10) 



Similarly the fcrmion two-point function is {tpiipj) = z i j . 
The vortex 4-point function can be computed from 



(T 12 T- 



■'SI/ 



1 r 
2 



^13^24 



214^23 



1/4 



- z 12 2: 23 2; 34 ;z 14 ' \ Z 13 Z 2A z l2 z 3i' 

(All) 

Taking the square root of this expression we get the 
known correlation funct ion of four cr's for which <j\ and a 2 
fuse to the identity!- 32 * 33 ! Now we use the conventions from 
the main part of the paper and take the limit AL — > 00. 
In this case only one of the terms in Eq. (All I survives 
and 



1 / z u z 23 



(Ti 2 T 3i ) — , 

AL^oo y/2 \Z X3 Z 2i Zi 2 Z 3A 



1/8 



We also have 



(^12> = 0, 



(A12) 



(A13) 



which is consistent with the notion that the tunneling of 
a vortex cannot create a fermion (or equivalently change 
the fermion parity). It is also straightforward to show 
that 



(tpiipjT^} = 0, 

AL— ^oo 



(A14) 



which means that a single vortex tunneling event is not 
enough to be able to transfer a fermion between the two 
edges. 



3. The six-point function 



(^jT 12 T 3i )(T[ 2 T 34 ) 

i) cos(4>j) COS 



2(cos( 



- j cos ( - 



4z i: j(z 12 z 34 ) 1/ 4 



j r/ ZgZ aZj2 Z 34 \ 1/2 + "I / £13*34 \ 



' Zl 3 Z 24 \ 1 /* 
•Z\ 4 Z 23 . 



Zi2Zi4ZjiZj 3 < 

f^Mi) 1/2 + (i^j)lf— Y'\ (A15) 

^Zi 2 Zi3ZjiZji/ i \Z-\zZ-)A/ > 



Dividing this with the square root of Eq. (All) the re- 
sult agrees with that of Ref. |331 We now take the limit 
of spatial separation AL — > 00, the only one term that 
remains is 



(A^T 12 T 3i ){T[ 2 T^) 
(z 13 z 24 ) 



(z i2 Z i 4Z j iZ j3 ) 1 / 2 



AL^oo ^ZilZ^Z^Zj^) 1 / 2 Zijiz^Z^ZuZ^) 1 / 4 



Combining this with Eq. (|A12| we find 

(ipiipjT 12 Tu) 



(z 13 z 24 r/* 



(Ai6) 



(Zi2Z i 4,ZjlZj 3 ) 1 / 2 



2 3 / 2 (z ll z i3 z j2 z :j 4) 1 /' 2 Zij (214^23) 1/2 

(A17) 

To get this result we have removed the phases associated 

with z 12 1 ^ 8 and Z34 . These phases are canceled when 
one makes sure that the tunneling term is described by a 
Hermitean term in the Hamiltonian. This is exactly the 
phase of (T 12 ) in Eq. ( |ATo| . 

Other orderings of the fcrmions and the tunneling 
terms are obtained by exchanging the indexes, for ex- 
ample 



{iiiT 12 xjjjT 3i ) 



(z 13 z 24 f/ s 



(z^Zi^ZijZj^ 1 / 2 



2 3 / 2 (z ll z i3 z 2j z :j 4) 1 / 2 (214223) 1/2 

(A18) 

The indexes on the z's should have the same order as 
they appear in in the original expression. This prescrip- 
tion was used in e.g. Ref. H5] and is equivalent to the 
Keldysh formalism for chiral bosons which is reviewed in 
e.g. Refs. [35] and [20j In the limit of spatial separation 
the last term gives a phase factor that depends on the 
order of the tunneling terms and the fermions according 
to 

p=ijl2M, 1234ij, 12ij34 
p = il234j (A19) 
p = il2j34, 12i34j. 




To calculate the contribution from a tunneling of two 
vortices we need the six-point function of two tp's and 
four cr's. This correlation function is a special case of the 
more general one that was first calculated in Ref. IMI with 
a similar method. To calculate the six-point function we 



APPENDIX B: EXCHANGE ALGEBRA 

An alternative formalism is provided by the exchange 
algebra of Ref. 1351 In this formalism the action of the spin 
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field is described by two types of operators a and b and 
their conjugates, a creates an excitation with dimension 
when acting on the vacuum, which is denoted by the 



shorthand a\0) — The conjugate <v interpolates 

in the opposite direction: a.T | — |q). Similarly b and 
frt interpolates between states of dimensions A and | 
according to b\jg) — ||) and b^\\) = |^). The exchange 
algebra is described by the following relations 



b\b 2 



V2 



a 2 a\ 
blh 



(Bl) 



which is Hcrmitean. 

Similarly we can represent the fermion field in terms 
of a's and 6's with coinciding arguments^ 



Mi, |0>-|A) 
>-|0) 



a[b[, 



(B6) 



Using Eqs. (B6) and (B3l together with the exchange 



algebra of Eqs. ( Bl I and (B2) it is straightforward to 



show that in all cases we have the following commutation 
relations 



a\a 2 
b x b\ 
b\a 2 



_ — «S12 



-«S12 -g g J Sl2 



tfct _ 



a\b 



(B2a) 
(B2b) 
(B2c) 
(B2d) 



The tunneling operator, e.g. Ti 2 , consists of a product of 
two er's in the identity channel, which we denote \<J\<j 2 \i. 
When acting on states with dimension or | this implies 
that we are allowed to use the representations 



[<J\o 2 \i 




(B3) 



Another important point is that the tunneling term 
should be represented by a Hermitean term in the Hamil- 
tonian. This can be achieved by explicitly adding the 
Hermitean conjugate in the definition of the tunneling 
term: 



T 12 cx [a x a 2 + a 2 (Tx] I = 



(1 



12 S)[cr 2 (7l]/ 



(l + e^f)[aiC7 2 ] / 



(B4) 



In the last step we used Eqs. (B3 1 and (B2 1. By adjusting 



the amplitude to conform with the result of the previous 



section [see discussion below Eq. ( A17 1] we define 

Tia = e-^ft^ao-i]! - e is ^[aia 2 ]i, (B5) 



(B7) 



which immediately implies the commutation relation be- 
tween a tunneling term and a fermion is 

Tuip 3 = si 3 s 23 ip 3 T 12 . (B8) 
With this very important relation we can always trans- 
form the correlation functions that we want to calculate 



(see Sec. VBl into one of two different forms: 



{■4>jT 12 T 34: Tpi}, 



(B9a) 
(B9b) 



Using the exchange algebra we can cluster decompose the 
last two expressions, in the limit of spatial separation we 
arc left with 



(ipjT 12 T 3 4ipi 



(ipjT 3 4T 12 ipi 



g-*Sl2 2 



(^■6 2 a 4 )(a t 1 6^ l ),(B10a) 



AL- 



AL- 



V2 



'3 l . 
,t|,t. 



V2 



-(^6 4 a 2 )(46 T 1 ^>(B10b) 



We have checked that the result of the formalism in this 
appendix gives identical results to those of the formalism 
in App. [X] Although the exchange algebra is derived at 
T = it also holds at finite temperatures. 
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